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We study the problem of a Dirac field in the background of an Aharonov-Bohm flux string. We 
exclude the origin by imposing spectral boundary conditions at a finite radius then shrinked to zero. 
Thus, we obtain a behaviour of the eigenfunctions which is compatible with the self-adjointness of 
the radial Hamiltonian and the invariance under integer translations of the reduced flux. 

After confining the theory to a finite region, we check the consistency with the index theorem, 
and evaluate its vacuum fermionic number and Casimir energy. 

I. INTRODUCTION 

The consideration of Bohm- Aharonov |jj scenarios is relevant to different physical problems, mainly 2+ 1-dimensional 
models in superconductivity and particle theory. The inclusion of the spin is an aspect to which much attention has 
been devoted. Many authors have studied this point in connection with the interaction of cosmic strings with matter 



We will treat the 3 + 1-dimensional problem of a Dirac field in the presence of a flux string. Due to the symmetry of 
the background field, the original four dimensional Dirac Hamiltonian can be written in a block-diagonal form, thus 
leading to two 2 + 1-dimensional problems. 

In this context, after developing the eigenfunctions in a convenient angular basis, the need to consider self-adjoint 
extensions of the radial Dirac Hamiltonian was recognized Q]^] . A one-parameter family of boundary conditions at 
the origin was then shown to arise, which amounts to establishing a relationship between components of the spinor 
rather than asking for the simultaneous finiteness of both of them. However, it was shown in |j] that only two values 
of the extension parameter correspond to the presence of a Dirac delta magnetic field at the origin. 

In references (2]]j]||, one of this possible boundary conditions was obtained, starting from a model in which the 
continuity of both components of the Dirac spinor is imposed at finite radius and the zero-radius limit is taken. 
As pointed out in ^,|l0] ^3|,when this boundary condition is imposed at the origin, the invariance under integer 
translations of the reduced magnetic flux (or, equivalently, large gauge transformations) is lost. 

In this paper, we adopt the view that the origin is an excluded point. The plane is thus a punctured one, which 
has the topology of a cylinder, and the aforementioned singular gauge transformations must constitute an invariance 
of the theory. It is in this spirit that, after devoting Sec. || to setting the problem, in Sec. Ill we impose on the Dirac 



fields spectral boundary conditions of the Atiyah- Patodi- Singer (APS) type ]14|-[l7[, taken as in reference |18| , |l9| 
at a finite radius. We show that, when the zero radius limit is taken one then obtains, at the origin, other among 
the boundary conditions compatible with self-adjointness, which present a minimal divergence behaviour pd] , ^2] , ^0| . 
They respect the above mentioned periodicity of the flux and are, at the same time, compatible with the presence of 
a Dirac delta magnetic field at the origin. 



From Sec. [V on, we confine the whole system within a punctured circle of radius R, imposing at the external 
boundary spectral conditions complementary to those satisfied at the position of the flux tube. We discuss the 
problem of zero modes and the compatibility with APS index theorem for manifolds with boundaries. 

In Sec. we determine the energy spectrum and study the fermionic number and Casimir energy of the confined 
system. This last calculation is performed in a Q function approach, following the lines of reference [p0[. 



Finally, Sec. VI contains some comments and conclusions. 



II. SETTING OF THE PROBLEM 

We study the Dirac equation for a massless particle in four dimensional Minkowski space. 
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(i 0-4)9 = (1) 
in the presence of a flux tube located at the origin, i.e., 

H = V AA = -S(r)e z (2) 
r 

where k = is the reduced flux. 

As the gauge potential is z— independent, equation (Q) can be decoupled into two uncoupled two-component equa- 
tions ||,by choosing: 

^° = (» °») ^ = {'°l) T * = ( "0" i°, ) ->> = (.",?) < 3 > 

In order to avoid singularities, we will consider that only r > tq is accessible, and take the limit ro — * 0, which is 
equivalent to having a punctured plane with the removed point corresponding to the string position. 
By taking 

A z = A r = A e = - , for r > r (4) 

r 

the Hamiltonian can be seen to be block-diagonal 

'-{"<? 1) o 

where the two-by-two blocks are given by 

/ ie* id (d r ±B)\ 

H± -\-ie± i6 (-d r ±B) ) (6) 

with 

B = --do - - (7) 
r r 

It should be noticed that these two "polarizations", which we will label with s = ±1, correspond to the two 
inequivalent choices for the gamma matrices in 2+1 dimensions 

From now on, we will be working with s — 1 (the case s = — 1 can be studied in a similar way, and explicit reference 
will be made to it whenever necessary). 

In this case, we can write: 

H+ = ( ° L q ^ , with L = -ie ie {-d r + B) Lt = ie -<s + B ) (8) 



and its eigenfunctions: 



9 E =( ^'2) , satisfy: %e = E Xe 



Now, the two components in 9e have different 9 dependence. In order to make this fact explicit, and to discuss 
boundary conditions at r = ro, we introduce JL SUl S|] 



ft V e^xiE {r,t 
and 

Li = -<9 r + B L[ = S r +B , (11) 

so that 
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L = -i e -^L x S-^ £t =i e —— iAe^y/i (12) 



and 



Li^ib = iE\iE 
l\xie = -iELpi E 

We expand tp\E and xis i n terms of eigenfunctions of B 



(13) 



(14) 



which are of the form: 

e n = e< n+ i) e , with A w (r) = - + ^ - , neZ (15) 



once the condition has been imposed that <ps and xe hi equation (^]) are single- valued in 6. 
Thus, we have 



XiE (r,«) = EL=oS->W«' W)fl 



(16) 



Replacing in (|13|) we obtain, for noninteger k (k = k + a, with k the integer part of k and a its fractionary part) 
Vie (r, 0) = V^E« (A. Jn-n (\E\r) + B n J K _„ (\E\r)) e< n+ i) e 



and 



Xie (r, 9) = -iJf v^E^-oo (A. ■/«+!-« (I^k) - J«-„-i (|£|r)) e <" + *) 9 ' '" 



xT , ,,, , Er=-oo(AJn-K(m+BnJ«-n(|^r))e- e 



En=-oo -i'-w (A. (|^|r) - B n J«-„-i (TO) e 

(Of course, for integer k, a linear combination of Bessel and Neumann functions must be taken). 
Finally, for s = — 1, the upper and lower components of ^e interchange, and E — > — _E. 



III. BOUNDARY CONDITIONS AT THE ORIGIN 

As is well known [|J|7]J^] , the radial Dirac Hamiltonian in the background of an Aharonov-Bohm gauge field requires 
a self-adjoint extension for the critical subspace n = k. In fact, imposing regularity of both components of the Dirac 
field at the origin is too strong a requirement, except for integer flux. Rather, one has to apply the theory of Von 
Neumann deficiency indices pl| , which leads to a one parameter family of allowed boundary conditions, characterized 

by§ 

i lim {MrT +1 g n (r) sin (j + f ) = M (Mr)"" f n (r) cos + |) (19) 

with v varying between —1 and [y — —a for s = X; u = a — 1 for s = —1) . Here, parametrizes the admissible 
self-adjoint extensions, and M, a mass parameter, is introduced for dimensional reasons. 

Which of these boundary conditions to impose depends on the physical situation under study. 

One possibility is to take a finite flux tube, ask for continuity of both components of the Dirac field at finite 
radius and then let this radius go to zero Thus, one of the possible self-adjoint extensions is obtained, which 

corresponds to O = fsgn (k) . 

As pointed out in j^|,[l0|,[l2| , this kind of procedure leads to a boundary condition that breaks the invariance under 
k — ► k + n (neZ). Now, this is a large gauge symmetry, which of course is singular when considering the whole plane, 
but is not so when the origin is removed or, equivalently, the plane has the topology of a cylinder. 

In order to preserve the aforementioned symmetry we propose, instead, to exclude the origin, by imposing spectral 
boundary conditions of the Atiyah-Patodi-Singer type 14 - 16) , as defined in |lq] , at a finite radius ro and letting 
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We consider the development in eq.(|i6|) and impose, at r = tq ■ 

fn (r ) = for A„ (r ) < , . 

g„M = for A„ (r ) > {M) 

for s = 1. 

As is well known, imposing this kind of boundary condition is equivalent to removing the boundary, by attaching 
a semi-infinite tube at its position and then extending the Dirac equation by a constant extension of the gauge field, 
while asking that zero modes be square integrable J^,^2| (except for A„ = 0, where a constant zero mode remains, 
with a nonzero lower component ). 

Then, going to eq.(|l7|) we have, after using the dominant behaviour of Bessel functions for small arguments : 



An 2{k+a-n) , 1 



n — k + — — a < 



D™ ^ 1 

-J- ~ V n-k + - 



(21) 



Now, we analyze two different situations: 



If a> | 



B n 



A n 

and the cigenfunctions in eq(|l8|) are of the form 



r D ->o for n < k 
r Q — »o for n > k + 1 



(22) 



* ( m V- r / J fe+Q _ n (|^r) e ^ \ ^ / J n _ fc _ Q (|£|r)e«* \ 

* ( ' } = A " I "i^^n-i ») e^ +1 > e J + j^i " I <J|1 J„ + r-,- Q (Wr) e««+* ) (23) 



n= — oo 



If a < \ 



A n 
B n 

An 



*y — >o for n < k — 1 
>r ^o for n > k 



(24) 



and the eigenfunctions are 

k-l 



**M) = ^ B„( (|^|r)e^J + 2- (,<Jf J B+1 _ fc _ a (|£?|r) e ^ J (25) 



n— — oo 



Notice that our procedure leads precisely to a self adjoint extension satisfying the condition of minimal irregularity 
(the radial functions diverge at r — ► at most as r _p , with p < h). It corresponds to the values of the parameter O : 

HI £-<i (26) 

As shown in O = ±^ are the only two possible values of the parameter which correspond to having a Dirac 
delta magnetic field at the origin. 

Moreover, this extension is compatible with periodicity in k. In fact, the dependence on k can be reduced to an 
overall phase factor in the eigenfunctions. 

As regards charge conjugation 

— > ai^E* ; « — > —k 

it is respected by the eigenfunctions, except for a = |. This is due to the already commented presence of a constant 
zero mode on the cylinder. 

Notice that, for the representation s = — 1 of 2 x 2 Dirac matrices, 
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** =u-*s>**>'j <27) 

and APS boundary conditions must be reversed, for A 7^ 0, as compared to (p0|), since the operator B changes into 
— B. For A = the lower component will be taken to be zero at ro which, as we will show later, allows for charge 
conjugation to be a symmetry of the whole model. Thus we take 

f n M = for A„ (r ) > , , 

5n (r ) = for A„(r )<0 

for s = — 1 

In this case, the resulting extension corresponds to 

This is an example of a physical application of APS boundary conditions, which are generally chosen due to their 
mathematical interest. We wish to stress that our conclusions concerning the behaviour at the origin hold true for 
massive Dirac fields. The procedure of imposing boundary conditions at a finite radius, then taken to zero, was 
considered for massive fields in [p3[ , where some vacuum quantum numbers were examined. 

It is worth pointing that, for integer k = k, our procedure leads (both for s = ±1) to the requirement of regularity 
of both components at the origin. In this case 



..Ml- £ Mfe^Mrt^l 



IV. THE THEORY IN A BOUNDED REGION. ZERO MODES AND APS INDEX THEOREM 



From now on, we will confine the Dirac fields inside a bounded region, by introducing a boundary at r = R, and 
imposing there boundary conditions of the APS type, complementary to the ones considered at r = ro. 
For s = 1, these boundary conditions are 

fn(R)=0 for X n (R)>0 . . 

g n (R) = for A n (R) < 

We start by studying the zero modes of our theory, which are of the form 

••<'■*- Ufe^*"*) (32) 



The conditions at r — rg imply 



A n — , for n < k + a - | 



B n = , for n>Ha 



Now, the boundary conditions at R imply 

A n — , for n > k + a - i , . 

B n = , for n < k + a - \ { ' 

Thus, no zero mode remains under this boundary condition, even without taking ro — > 0. This is in agreement with 
the APS index theorem [^4|-^6|. In fact, according to such theorem 

n + -n-=A+b(r )+b{R) (35) 

where n+(n_) is the number of chirality positive (negative) zero energy solutions, A is the anomaly, or bulk contri- 
bution, and b are the surface contributions coming from both boundaries |l7jjil| 
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b(R) = -(h R -r)(R)) 



b{ra) = g 0?M ~ h ro ) 



(36) 



with 



V(r)= ^2 s S nA «M \ x n(r)\~ 

A„(r)^0 



(37) 



s=0 



the spectral asymmetry of the boundary operator B and h r is the dimension of its kernel. 
In our case : 



b(ro) 




HR) 




a > 



a < 



(38) 



Thus, the boundary contributions cancel. As regards the volume part, it also vanishes for the gauge field configu- 
ration under study, and we have n + — n_ = 0, which is consistent with our analysis following eq. (|34|) . 
For s = — 1 both boundary contributions interchange, and identical conclusions hold regarding the index. 



V. ENERGY SPECTRUM. FERMIONIC NUMBER AND CASIMIR ENERGY 

The energy spectrum can be determined by imposing (for s = 1) the boundary conditions ( |3l"| ) at r = R on the 
eigenfunctions of eq.(p3|) if a > h, or on the eigenfunctions in eq. J25|) , if a < 4. When doing so, we obtain : 

± J "^°'' , n=l, ...oo Z = l...oo r ^ 1 

±a 

R 



= ' fora>- (39) 



and 



I ± J " ,n = 0, . . .oo Z = 1. ..oo r 1 

E ",i = 1 , j^ti n /i fora<- 40 

± t, n = Q .oo t = l...oo 2 



where .?'„,; is the Z-th positive root of Jj/. The same spectrum results for s = — 1. 

For both s values, the energy spectrum is symmetric with respect to zero. This fact, together with the absence of 
zero modes results in a null vacuum expectation value for the fermionic charge pjj ] : 

<N>+=-±(n + -n_) = (41) 
For the same reasons < N >_= 0, so the total fermionic number of the theory is null. 

It is interesting to note that the origin contributes to the fermionic number with (see the discussion in p3fl) 



< N > ro ,±= { ±\ a=\ (42) 




which coincides, for each s value, with the result presented for the whole punctured plane in |10|,|12|, where opposite 
signs of the mass correspond to our opposite signs of s (see also [|6| for related work) , except for a = h. This last 
fact is associated with charge conjugation non invariance in each subspace. However, the sum of both contributions 
cancels for all a. 

We go now to the evaluation of the Casimir energy, which is formally given by 

(43) 



\E>0 E<0 ! \e>Q 



where the symmetry of the energy spectrum for our problem has been used. Of course, a regularization method must 
be introduced in order to give sense to this divergent sum. In the framework of the C-regularization p7|p8| , 



G 



where the parameter fi was introduced for dimensional reasons. 

Here, it is useful to define the so-called partial zeta function as in references po|,P9|.pO| 



(44) 



c(*) = X>-*,i)- 



1=1 



So, for the problem at hand we have, from eqs. (|39|) , ([io|) , (|44|) and ME 



(45) 



with 



Sn=-oo C|n-a| + Cq-1 
oo C|n— a| + C— a 



2=— 1 



for 
for 



a > 
a < 



(46) 



(47) 



Notice that, as a consequence of the invariance properties of the imposed boundary conditions, the Casimir energy 
is both periodic in k and invariant under a — > 1 — a, as well as continuous at integer values of k. 

For any value of K, it is the sum of the energy corresponding to a scalar field in the presence of a flux string and 
subject to Dirichlet boundary conditions (studied in ref. |p0| ), plus a partial zeta coming, for fractionary k, from the 
presence of an eigenfunction which is singular at the origin or, for integer k, from the duplication of Jo- 

Both contributions can be studied following the methods employed in P0| and developed by the same authors in 
previous work p9|p0[ . So, here we won't go into the details of such calculation. 

The scalar field contribution presents a pole at z — — 1, with an a-independent residue. So, in this case, a 
renormalized Casimir energy can be defined as 



E c (a) - E C (0) 



(48) 



Following PQl , this quantity can be obtained through a numerical calculation. (Notice there is a minor error in 
eq.(3.6) of that reference. In fact, a factor of tt is missing in the fifth term on the right hand side of such equation). 

Now, going back to the Dirac field, the partial zeta appearing in (^6|),(p7|) can also be evaluated with the methods 
of pot]. However, it presents a pole, whose residue is a-dependent : 



1 

8^ 



l-4z/ 



1 



1 



(49) 



Notice that such residue vanishes only for v — ±| (a = \ in the case of interest). 

The finite part of the Casimir energy is plotted in fig.(p|) as a function of a for /iR =1. Of course, an absolute 
meaning cannot be assigned to it, due to the presence of the pole and the consequent need to introduce a-dependent 
counterterms. 



VI. COMMENTS AND CONCLUSIONS 



In summary, we have studied the problem of Dirac fields in the presence of a Bohm-Aharonov background which 
can, by virtue of the symmetry properties of the external field, be separated into two 2+1-dimensional ones (labelled 
through the paper by s — ±1). 

We have adopted the viewpoint that the plane is a punctured one, so that the model must turn out to be invariant 
under integer translations of the reduced flux, which are nothing but "singular" gauge transformations. In this spirit 
we have shown that, by imposing spectral boundary conditions of the APS type at a finite radius and then letting this 
radius to zero, one among the family of allowed boundary conditions at the origin is obtained for the two-component 
spinors. Under this boundary condition, the dependence on the integer part of the flux can be extracted from the 
eigenfunctions of the Hamiltonian as an overall phase. As a consequence, the aforementioned symmetry is preserved. 

As regards the invariance under charge conjugation, it is broken for each s value by the boundary conditions at the 
origin, for the fractionary part of the flux equal to \. This is due to the fact that a constant zero mode of definite 
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handedness is allowed for on the cylinder when applying the APS boundary conditions. However, this symmetry is 
preserved when opposite handed conditions are imposed on both s-subspaces. 

It is worth stressing that this approach to the problem of self-adjointness at the origin is, to our knowledge, the 
first proposal of a physical application of APS boundary conditions in this context. 

We have then confined the theory to a finite region, i.e. a punctured circle, and imposed adequate APS boundary 
conditions at the exterior boundary, with the aim of checking consistency with the APS index theorem for manifolds 
with boundaries, which we have shown to hold. 

Finally, we have studied the vacuum fermionic number and analized, through zeta function regularization techniques, 
the Casimir energy of the model. Since this last turns out to present a divergence depending on the fractionary part 
of the flux, no absolute significance can be assigned to its finite part, which is thus subject to renormalizations. 

However, notice that, contrary to the treatment of the origin, APS boundary conditions were imposed at the exterior 
boundary for merely formal reasons (the existence of an index theorem for this case) . The vacuum energy under local 
(bag-like) external boundary conditions is at present under study. 
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FIG. 1. Finite part of the adimensionalized Casimir Energy 
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FIGURE CAPTION 

FIG. 1 - Finite part of the adimensionalized Casimir Energy 
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